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Relativistic, spherically symmetric configurations consisting of a gravitating magnetized 
anisotropic fluid are studied. For such configurations, we obtain static equilibrium solutions with an 
axisymmetric, poloidal magnetic field produced by toroidal electric currents. The presence of such 
a field results in small deviations of the shape of the configuration from spherical symmetry. This 
in turn leads to the modification of an equation for the current and correspondingly to changes in 
the structure of the internal magnetic field for the systems supported by the anisotropic fluid, in 
contrast to the case of an isotropic fluid, where such deviations do not affect the magnetic field. 
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I. INTRODUCTION 


Magnetic fields are widely believed to play an important role in modeling compact, strongly gravitating objects - 
neutron stars. The measured surface magnetic fields of such objects are in the range of 10 11 — 10 13 G (for the so-called 
“classical pulsars”) and can reach values of the order of 10 15 G for magnetars. Considerable research has been devoted 
to stellar magnetic fields (see, e.g., the book [I] on stellar magnetism in general and some works on particular 

questions of the structure of strongly magnetized configurations, and also references therein). The studies show that 
strong magnetic fields affect considerably the interior structure of neutron stars and determine their evolution in time 
(the magnetic dipole radiation, the deformation due to the magnetic stress, etc.). 

In modeling magnetized neutron stars, neutron matter is usually described as being an isotropic perfect fluid, i.e., 
a fluid obeying Pascal’s law. However, such a description cannot be considered as completely satisfactory, since at 
extreme densities of matter and magnetic fields typical for the interior regions of neutron stars, it is already necessary 
to take into account the interaction of matter with ultra-strong magnetic fields. The presence of such fields may lead 
to the appearance of unequal principal stresses in the neutron fluid (see Ref. [|J and references therein); i.e., the fluid 
will become anisotropic. (See, however, Ref. [b| for further discussion of the question of whether the hydrodynamic 
pressure is isotropic or anisotropic in the presence of magnetic fields.) Among the other possible reasons for the 
appearance of the anisotropy in superdense matter might be the nuclear interactions [f|, pion condensation |7|, some 
kinds of phase transitions [U, or viscosity [S|. 

Regardless of the specific nature of the fluid anisotropy, its presence may result in considerable changes in the char¬ 
acteristics of relativistic stars, as demonstrated, for example, in Refs. EMI- But apart from that, when considering 
magnetized matter, one might also expect the influence of the anisotropy on the structure of the magnetic field. The 
aim of the present paper is to get some insight into this question. 

We emphasize that our goal is not in itself to model more or less realistic distributions of the magnetic field. We 
just want to demonstrate some effects associated with the presence of the fluid anisotropy. To do this, we employ a 
simplified model of a poloidal magnetic field, i.e., the field lying in the meridional planes. Such a field is known to 
be unstable 13 15|, and only configurations consisting of poloidal and toroidal fields with comparable strengths can 


16| . Nevertheless, the model with a purely poloidal field, on the one hand, reflects the main qualitative 


be stable m 

characteristics of more realistic magnetic field structures, and, on the other hand, within this model it is quite easy 
to carry out calculations and to keep track of changes coming from the fluid anisotropy. 

The presence of such a dipole field will inevitably lead to deformations of the configuration, so that it is no longer 
spherically symmetric but possesses instead axial symmetry. However, since the energy density of the magnetic field 
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discussed here is much smaller than the energy density of the neutron matter, we follow Ref. E2 and consider these 
deformations as second-order corrections to the background spherically symmetric configuration. Using this fact, 
we make a perturbative expansion of the field equations. In this case the equations for fluid perturbations are not 
independent but subject to an integrability condition, which in turn gives an equation for the current. Its solution 
depends considerably on the degree of fluid anisotropy, and this ultimately causes changes in the structure of the 
magnetic field of the configurations under consideration. Note that this approach was used by us earlier in Ref. [xl] 
to describe magnetic fields in mixed systems consisting of a wormhole filled by a strongly magnetized isotropic or 
anisotropic neutron fluid. 

The paper is organized as follows: In Sec. [II] we present the statement of the problem and derive the corresponding 
general-relativistic equations for the systems under consideration. In Sec. IHIl we numerically solve these equations for 
isotropic and anisotropic cases. Comparing the results, we demonstrate the influence of the fluid anisotropy on the 
structure and strength of the magnetic fields. Finally, in Sec.[|V]we summarize the results obtained. 


II. STATEMENT OF THE PROBLEM AND EQUATIONS 

Here we consider gravitating configurations consisting of a strongly magnetized anisotropic fluid. Our purpose is 
to examine how the presence of the fluid anisotropy affects the structure of the magnetic field. In doing so, we make 
use of the following simplifying assumptions 0: ' 

• The magnetic field is modeled in the form of an axisymmetric, poloidal magnetic held produced by toroidal 
electric currents. Its presence in the system will in general result in a deviation of the shape of the configuration 
from spherical symmetry. However, for values of the magnetic held strength of the order of 10 12 — 10 15 G, which 
will be discussed below, deviations from the spherical shape are negligible, since the energy of the magnetic 
held is much smaller than the gravitational energy {*191] . This permits us to neglect in the zeroth approximation 
the deformations of the conhguration associated with the magnetic held and to consider such deformations as a 
second-order effect. 

• We do not take into account rotational deformations and consider only static equilibrium conhgurations. 

• The neutron matter is modeled by the simplest polytropic equation of state. 

• The interior of a star is assumed to be a perfectly conducting medium free of electric charges and helds. 

Consistent with this, in Sec. Ill Al we write down the background equations for a spherically symmetric case, and in 
Sec. Ill Bl we derive an equation for the magnetic held. Finally, in Sec. Ill Cl we perturb the background solutions and 
employ the obtained perturbed equations to derive an equation for the current. 


A. Background equations 

To derive the Einstein equations and the Tolman-Oppenheimer-Voikoff equation for the huid, we choose the spher¬ 
ically symmetric line element in the form 

ds 2 = e v (dx 0 ) 2 - e x dr 2 - r 2 (dQ 2 + sin 2 0 def 2 ), (1) 

where v and A are functions of the radial coordinate r only, and x° = ct is the time coordinate. We refer to these 
equations as the background equations. 

As a matter source in these background equations, we take an anisotropic huid, i.e., the huid for which the radial, p r , 
and tangential, p t , components of pressure are not equal to each other. For such a huid, the energy-momentum tensor 
can be taken in the form (see, e.g., Ref. fl3j l 

T Ufl) = ( £ + Pt) u ^ u v “ S uPt + (Pr - Pt) S^S U , (2) 

where e is the huid energy density. The radial unit vector s M is dehned as = (0, e -A / 2 , 0,0), with = — 1 

and = 0. The energy-momentum tensor then contains the following nonzero diagonal components: = 

(e,-Pr,-Pt,-Pt)- 
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Using the metric (HD and the energy-momentum 
Einstein equations, 


tensor m, one can obtain the (*) and (£) components of the 
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where G is the gravitational constant and the prime denotes differentiation with respect to r. In turn, the p = r 
component of the conservation law, T= 0, yields the generalized Tolman-Oppenheimer-Volkoff equation for the 
anisotropic matter: 


dp r nl , M + 4nr 3 p r / c 2 2 

— = -G (e + p r ) —— - +-{Pt~ Pr) ■ 

dr r(c z r — 2GM) r 

In obtaining this equation, we have used Eq. © and introduced a new function M(r), defined as 

z -x._ 1 2 GM(r) ' 


Substituting this expression into Eq. m, we also have 


dM 47t 2 
dr c 2 


( 5 ) 


( 6 ) 


For a complete description of the background configuration, Eqs. © and © must be supplemented by an equation 
of state (EOS) for the fluid. Here we restrict ourselves to a simple barotropic EOS where the pressure is a function of 
the mass density p. In this case, there are two possibilities of specifying the EOS. First, one may take two different 
EOSs for the radial and the tangential components of pressure, p r = p r {p) and p t = Pt(p)- Second, it is possible to 
restrict oneself to one EOS, say, p r = p r {p ), but, in addition to this, one may then assign the function A = p t — p r 
which appears in Eq. ©• This function is called the anisotropy factor j20|. 

Here we employ the second possibility, using for the radial pressure the simplest polytropic EOS in the form 


p r = Kp 1 , e = pc 2 , 


( 7 ) 


where K , 7 are constants whose value depends on the properties of the fluid under consideration. Such an EOS was 
used in Ref. 0 in modeling general relativistic isotropic fluid spheres and in Refs. 0, |lj|] in describing objects with 
an anisotropic pressure. 

Following Ref. [TH , here we take the anisotropy factor A in the following form: 

A = pt-p r = ap r p , ( 8 ) 


where a is the anisotropy parameter, and the function 

2 GM{r) 


is called the compactness. 

The choice © has the following attractive features 0 : (i) since at r —> 0 the compactness p ~ r 2 , the anisotropy 
factor vanishes at the center (where the fluid becomes isotropic), and this ensures the regularity of the right-hand side 
of Eq. © (for other possibilities to obtain regular solutions without imposing the requirement for the anisotropy to 
vanish at the center, see Ref. 0 ); (ii) the anisotropy factor thus defined makes itself felt only for what are essentially 
relativistic configurations, for which p ~ G>{\). This agrees with generally accepted folklore, according to which the 
fluid anisotropy may play an important role only at high densities of matter [10-Il2l|. 

Note that in the absence of p on the right-hand side of © the ratio of pressures pt/p r would be constant along the 
radius, and the problem would become analogous to the case of Refs. [ 10, 1131 • However, such an ansatz for the fluid 
anisotropy seems to be less realistic than that of ©. Another ansatz for the anisotropy can be found in Refs. 0 , 0 . 
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B. Equation for the magnetic field 


As an ansatz for the magnetic field, choose an axisymmetric, poloidal magnetic field, which is created by a 4-current 
j M =(o,o,o,j» Ea. For such a current, the electromagnetic 4-potential A^ has only a ^-component A M = (0,0, 0, A^). 
In this case, taking into account the nonvanishing components of the electromagnetic field tensor F r( f, = dA^/dr and 
Fq ^ = dA^/dO, the Maxwell equations give in the metric (HJ the following elliptic equation: 




1 d 2 A (j) 1 

r 2 ,902 r 2 CO U d Q 


= - c u. 


(9) 


The solution of this equation is sought as an expansion of the potential A$ and the current j<f, as follows da Hi; 


A# = ^a;(r) sin 0 


1=1 

OO 


dPi( cos0) 
d& 5 


U = sin 0 


i=i 


dPi (cos 0 ) 
d@ ' 


( 10 ) 

( 11 ) 


where Pi is the Legendre polynomial of degree l. Substituting these expressions into Eq. ED, we have 

e-V + \ («/ - A') e~ x a[ - = -~j t . (12) 

2 r z c 

The solution of this equation is sought for a given current ji, an equation for which will be derived in the next 
subsection. Since here we consider only a dipole magnetic field for which l = 1, for convenience, we drop the 
subscript 1 on a and j hereafter. 


C. Integrability condition 

For the magnetic field under consideration, the current j is not arbitrary but must satisfy an integrability condi¬ 
tion @,113- To obtain this condition, we make use of the fact that the magnetic field induces only small deviations 
in the shape of the background spherically symmetric configuration. To describe these deviations, let us employ the 
approach adopted in Ref. da and expand the metric in multipoles around the spherically symmetric spacetime. Then 
the deformations of the metric and the fluid are regarded as second-order perturbations, whereas the electromagnetic 
field potential and the current are regarded as first-order perturbations. In this case the corresponding metric can be 
taken in the form 

( O p A(r) 'l 

ds 2 = e{1 + 2 [ho(r) + h 2 (r)P 2 {cos 0 )]} ( dx °) 2 — e A ^ s 1 + —-[m-o( r ) + m 2 {r)P 2 (cos 0 )] > dr 2 

—r 2 [1 + 2fc 2 (r)P 2 (cos0)] (c?0 2 + sin 2 0d^ 2 ), (13) 

where ho, h 2 , mo, m 2 , and k 2 are the second-order corrections of the metric, and P 2 denotes the Legendre polynomial 
of order 2 . 

The total energy-momentum tensor for the system under consideration is 

T$ = (e + p t ) u»u v - Sipt + (p r - Pt) s^s„ - + -^F aft F^. (14) 

Using the relation between the pressures ([5p. we can eliminate the tangential pressure; i.e., henceforth we will work 
only with the radial pressure. 

Expand now the fluid energy density and the pressure in the form 

e(r, 0 ) = £ 0 + ^7 {Spo + SP 2 P 2 ), (15) 

Po 

p(r, 0 ) =p 0 +Sp 0 + Sp 2 P 2 , (16) 

where the background solutions Sq,Po and the perturbations Spo, Sp 2 depend on r only. (For convenience, we hereafter 
drop the subscript “r” on p.) Substituting these expressions into the conservation law, Tjf = 0, and using the metric 
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(usd , we obtain the following v = r and v = 0 components: 


5p' 2 = 2app 0 k 2 - (e 0 + p 0 ) h 2 - - 


u' ( 1 + ^ ^ 

Po 


A 2 «' 3 
<>P 2 - o ~2 
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{l + ap)dp 2 = - [{1 + ap)p 0 + £ 0 \ h 2 - -- o~2~- 

1 — p r ore 

The function v' entering into Eq. El) can be found from the background equations in the form 

p , = 2 p' 0 f 2appo _ A 

£o+Po\ r p' 0 ) ' 

The integrability condition for Eqs. El and (fl 8 l) gives the following equation for the current: 

j' + Fj + N = 0 

with 


F = -- 
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(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 


Note that all terms in Eq. © are of first order. Indeed, j itself is of first order, and the function F is of zeroth order. 
Finally, the function N is also of first order, since the potential a, which is of first order, enters in the denominator, 
and the numerator contains the second-order terms. 

In the anisotropic case (a ^ 0), the function N is nonzero. To evaluate it, one needs to find the functions h 2 and 
k 2 appearing in (12211 . which are determined by the following set of equations deduced from the Einstein equations: 


,, / r a , 1 [3 2 1 + rz/ 
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(23) 

(24) 
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Notice here that in deriving the expression ( 1221 ) . the function m 2 , which appears in (1181) . has been eliminated by 
using Eq. (EH). 

For the isotropic case (a = 0), Eq. (E0l) can be integrated analytically to give 


j = c 0 r 2 (£ 0 +po), (27) 

where Co is an integration constant. This expression coincides with the one obtained in Ref. E3 [see their Eq. (25)]. 

It should be emphasized that the small deviations in the shape of the configuration described by the second-order 
corrections to the metric h 2 , k 2 , and m 2 affect the current, and correspondingly the structure of the magnetic field, 
only if the fluid anisotropy is present. The question of how this manifests itself will be discussed in the next section. 


III. NUMERICAL RESULTS 

In this section we numerically integrate the obtained equations for isotropic ( a = 0) and anisotropic (a ^ 0) cases. 
Our goal is to exhibit differences in distribution of the magnetic field for these two cases. 
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A. Background equations 


For numerical calculations, it is convenient to rewrite Eqs. © and © in terms of dimensionless variables. To do 
this, we employ the usual reparametrization of the fluid density (23} , 


= p r 6 n 


p = p< 


(28) 


where p c is the density of the fluid at the center of the configuration, and the constant n, called the polytropic index, 
is related to 7 from the EOS © via n = 1/(7 — 1). Introducing the dimensionless variables 


f = j. «(0 


M(r) 
Airp c L 3 ’ 


where 


L = 


I a(n + l)c 2 
AnGp c 


we have from Eqs. © and © 


, , x dO 
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fe n . 


—a(n + 1)(1 + ad) 


v + a0?9 n+1 

e (1 - p) 


9 P® 
2 aa—, 


(29) 


(30) 

(31) 


Here a = Kp)/ n /c 2 = p c /(p c c 2 ) is a constant related to the pressure p c of the fluid at the center, and p = 2a(n + l)u/£. 
The parameter a , also called the relativity parameter [21], measures the deviation of a configuration from nonrela- 
tivistic systems for which a <C 1. 

These equations are to be solved for given n and K subject to the boundary conditions given in the neighborhood 
of the center by the following expansions: 

9 ~ 1 + 2 ^ 2 ? 2 ; v ~ where d 2 = — (a + l)(a + 1/3) + -aa. (32) 

Notice that in the case of an isotropic fluid (a = 0), 62 is always negative, and correspondingly the fluid has a 
maximum density at the center. For the anisotropic fluid (a ^ 0), a situation may occur where 6 2 > 0; i.e., the fluid 
will have a local minimum of the density at the center, and a maximum will be located somewhere between the center 
and the edge of the star. An example of such a configuration can be found in Ref. [T3 . 

Using these boundary conditions, we numerically solved Eqs. m and started the solutions near the origin 
(i.e., near £ « 0), and solved out to a point £ = £&, where the function 9 became zero (the boundary of the fluid where 
P(€b) = 0 ). 

Note here that the radial coordinate r from the metric © describes the areal radius of a sphere with area Anr 2 . 
Another physically relevant radial coordinate is given by the coordinate £ associated with the proper radius, which 
represents the true distance from the center. It is defined as follows: 


£ = 


fe^= f 

J 0 -'o 


da 


V 1 - p 


(33) 


Then, the proper radius of the fluid R is obtained in dimensional variables as R = £,bL- 


B. Magnetic field equations 


We rewrite now the perturbative equations (fl2l) and (l2i)l) for the electromagnetic potential and the current through 
appropriate dimensionless variables: 
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j' + Fj + N = 0, 
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where F, N are the dimensionless expressions obtained from F and N , Eqs. (EH) and E2l) . 

The boundary conditions for these equations are given by an expansion in the neighborhood of the center, 


a 


1 

2 





(37) 


where a c ,j c are free parameters. 

In the isotropic case N = 0, and Eq. (l36l) has the solution in the form of (1271) . 

When an anisotropy is present, Eqs. © and (1361) are solved together with the Einstein equations (1231) - (1261) for the 
second-order perturbations. To obtain regular solutions, we choose the boundary conditions for them in the form 

h 2 ~ \h 2 ct 2 , k 2 ~ \k 2c ^ 2 . (38) 

In turn, the boundary condition for the function m 2 can be found from Eq. (1251) . Using these boundary conditions in 
Eq. EH) or in Eq. EH) , one can find the following constraint on the expansion parameters: 


h 2c + k 2c = -a 


2 

C ‘ 


Thus, we have three free parameters, a c ,j c , and h 2c or k 2c , which are chosen in such a way as to (i) get a required 
value of the surface magnetic field, (ii) provide regularity of the perturbed solutions at the center, £ = 0, and (iii) obtain 
asymptotically decaying solutions for £ —> oo. 

In particular, it is necessary to derive the known external solution pa, according to which beyond the fluid the 
electromagnetic field potential decays as 


a 


-r 


In (l-n) + n+ -i? 


In this solution, /i corresponds to the external vacuum solution for the background configuration with the mass 
concentrated inside the radius £{,. 

Finally, by solving the equation for the electromagnetic field potential, we can find the strength of the magnetic 
field, which is given by the following tetrad components (i.e., the components measured by a locally inertial observer): 


Bf — Fq£ — ca 


2 p c 


er(n + 1) 


■ COS0, 


B t \ = F-i = —c 

0 rq> 


Pc 


2er(n + 1) £ 


v 71 - P -/ • 0 
a sm0. 


(39) 


C. Structure of the magnetic field 

Consider now magnetic field distributions for the configurations with different values of the anisotropy parameter a. 
In Ref. [l2j], finite-size, regular solutions for compact relativistic configurations with the EOS 0 and —2 < a < 4 
have been found. For such systems, the mass-radius curve has a turning point corresponding to the maximum allowed 
mass at which the system is still stable. It was shown that as a decreases, the maximum mass decreases as well. In 
particular, when the polytropic index n— 1 (7 = 2) and I\ = 6.67 x 10 4 cm 5 g -1 s~ 2 (in our units), the maximum 
mass M max ss 1M 0 for the system with a = —2 and M max « 2.14 Mq for the case of a = 4. 

Here we use the configurations of Ref. 0 as background systems in which the magnetic field is added. In order 
to show explicitly the influence of the fluid anisotropy on the structure and strength of the magnetic field inside the 
star, we will proceed in the following way: 

(i) Following Ref. Q, we fix the number of baryons contained in each of the configurations under consideration 
and compare the configuration supported by an isotropic fluid with the anisotropic systems. This approach implies 
that systems with the same baryon number (or equivalently, the same baryon mass Mb) may have different physical 
characteristics (total masses, sizes, an internal structure, etc.), depending both on the type of EOS and the presence 
of ultra-strong magnetic fields Q and on the presence of fluid anisotropy. As will be shown below, the presence of 
the latter will in turn lead to considerable changes in the structure of the magnetic field. 

(ii) All the systems in question are modeled by the EOS 0 with some fixed polytropic parameters. This allows us 
to keep track of changes in the structure of the magnetic field coming from just the fluid anisotropy, but not those 
which depend on the form of EOS. 

(iii) Depending on the value of the anisotropy parameter a, the total mass M of the configurations will change. 
(Note that the total mass is understood here to be only the mass of the background configuration without taking 









B-JB, 




j,x 10 5 



FIG. 1: The tetrad components B? and Bq of the magnetic field (in units of the surface strength of the magnetic field B s at 
the pole) evaluated on the symmetry axis (0 = 0) and in the equatorial plane (0 = tt/ 2), respectively, and the dimensionless 
current j are shown as functions of the relative radius £/£i>. The numbers near the curves denote the values of the anisotropy 
parameter a. The thin vertical lines correspond to the boundary of the fluid. All configurations have the same baryon mass, 
Mb ~ 1.23JV/q (for other characteristics, see Tabic [I]). For a = —1.7, the central values of Bf/B s and Bq/B s are of the order 
of 44 (modulus). 



FIG. 2: Magnetic field lines for the configurations of Fig. [T] The plots are made in a meridional plane </> = const, spanned by 
the coordinates x = rsin0 and z = rcos0. The circles denote the boundary of the neutron fluid, possessing the radius R. 
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FIG. 3: Energy density = £ of the neutron fluid (in units of p c c 2 ) and the ratio Pt/Pr of the tangential pressure to the 
radial pressure for the configurations of Fig. [T] For the left panel, the values of the anisotropy parameter a are successively 
decreased from a = 4 (for the top curve) to a = —1.7 (for the bottom curve). 


TABLE I: Characteristics of the configurations exhibited in Fig. [T] Here the central mass density p c of the neutron fluid (in 
units of 10 14 g cm -3 ), the relativity parameter a, the total mass M (in solar mass units), and the proper radius of the fluid R, 
as given by Eq. ([33]). (in kilometers) are shown. 


a 

pc 

a 

M/Mq 

R 

- 1.7 

48.00 

0.36 

1.02 

9.44 

-1.0 

27.60 

0.20 

1.05 

10.29 

0.0 

17.50 

0.13 

1.07 

11.12 

1.0 

12.78 

0.10 

1.08 

11.77 

2.0 

10.02 

0.07 

1.09 

12.31 

4.0 

6.92 

0.05 

1.10 

13.22 


account of the mass coming from the magnetic field that is negligibly small compared with the mass of the fluid.) In 
order to obtain the required fixed baryon mass Mb for the systems with different a, we will appropriately change p c . 
In turn, the free parameters of the magnetic field a c ,j c will be chosen in such a way as to provide the required surface 
magnetic field B s (say, at the pole) for all configurations, and also to obtain asymptotically vanishing solutions for 
the electromagnetic potential. 

Proceeding in this manner, we obtained the results shown in Figs. [Hand[2] All calculations were carried out using 
the polytropic parameters n = 1 (7 = 2) and K = 6.67 x 10 4 cm 5 g - 1 s -2 . Since for these parameters the system 
with a = —2 possesses the lowest maximum total mass M max ss 1M 0 , we took it as the reference mass. Then, we 
calculated the baryon mass Mb of this system, fixed this value and then calculated the corresponding total masses for 
the systems with different a > — 2 but all having this fixed value of Mb . The results of the calculations are listed in 
Table [J] Note that all the configurations thus obtained are energetically stable (cf. Ref. 0 )- 

In Fig. [1] the distribution of the component Bf on the symmetry axis (0 = 0), the distribution of the component 
Bq in the equatorial plane (0 = 7 r/ 2 ), and the current j are shown. It is seen that the magnitude and the distribution 
of the internal magnetic field change considerably for the systems with different a. In the case of a > 0 (a < 0), 
the value of the component Bf. is always less (greater) than that of the isotropic case. In turn, the magnitude of the 
component Bq for the anisotropic fluid may be either greater or less (modulus) than the one for the isotropic case, 
depending on the point where a comparison is performed. It is interesting to note that for positive values of a the 
profiles of both components Bf and Bq are flattened over the radius up to the point ~ i?/2. On the other hand, for 
large negative values of a, there is a rapid growth of the central values of Bf and Bq, which, for a = —1.7, differ in 
order of magnitude from those of the isotropic case. 
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The distribution of the current also depends considerably on the value of a, see Fig. [l] For example, if in the 
isotropic fluid the current always has a fixed sign [cf. Eq. (1271) ]. then in the anisotropic case it may have variable sign 
both for positive and for negative a. 

The magnetic field lines are shown in Fig. [2j where the spatial coordinates are given in units of the radius of the 
fluid, R. The behavior of the field lines is qualitatively similar for different values of a. 

Note that the perturbed equations (USD, (HDD-® are invariant under the transformation a,j —> (3a,/3j and 
h 2 ,k 2 ,m 2 —> (3 2 li2, (3 2 k2, /3 2 rri2, f3 being any constant. Correspondingly, the components of the strength of the mag¬ 
netic field given by Eq. (l39l) transform as B^,Bq —> (3Bf,(3BQ. Then, if one simultaneously replaces B s by (3B S , the 
plots shown in Fig. |T]are unchanged for any value of the field, and the dimensional values (in gauss) of the strength 
of the magnetic field are obtained on multiplying these plots by the required surface value B s . Of course, this holds 
true only in the approximation used here when one can neglect the influence of the magnetic field on the background 
solutions. In particular, the obtained plots are applicable both to the “classical pulsars” (for which B s ~ 10 12 G) and 
to magnetars (for which B s ~ 10 15 G). 

It may also be noted that for the systems with negative a, the relativity parameter a is greater than that for the 
case of positive a (see Table U. That is, the smaller a, the more relativistic matter is necessary in order to have 
equilibrium configurations with a fixed number of baryons of the type we discuss here. In turn, the growth of the 
relativity parameter a is accompanied by the following effects: (i) a greater concentration of matter toward the center, 
as demonstrated in the left panel of Fig. [3] where the distributions of the fluid energy density are shown; (ii) for 
strongly relativistic systems with large negative a there exists the above-mentioned rapid growth of the strength of 
the magnetic field in the internal regions of the configurations, see Fig. [TJ (iii) for negative a , a more rapid increase 
in the difference between the tangential and radial pressures takes place (as one can see, e.g., from a comparison of 
the curves for a = 1 and a = — 1 shown in the right panel of Fig. 0. 


IV. CONCLUSION 

We have studied equilibrium, gravitating configurations consisting of a strongly magnetized fluid with an anisotropic 
pressure. The fluid is described by the simplest polytropic EOS 0, and the anisotropy is modeled by the expres¬ 
sion (|g|). which takes into account both the local properties of the matter (through the pressure) and the quasilocal 
properties of the configuration (through the compactness). Our goal was to clarify the question of how the presence 
of such an anisotropy influences the structure of the magnetic field, which was modeled here in the form of a poloidal 
field produced by toroidal electric currents. 

For this purpose, we compared configurations with the same number of baryons but with different degrees of fluid 
anisotropy. Having fixed the polytropic parameters, we kept track of changes in the internal structure of the magnetic 
field coming from just the fluid anisotropy (for details, see Sec. lIII Cl) . These changes were then revealed by comparing 
the configurations with isotropic and anisotropic fluids. 

In summary, the main results of the studies are as follows: 

(i) For an anisotropic fluid, an equation for the current is modified in such a way that it begins to “feel” the 
presence of small deformations of a star associated with the presence of a dipole magnetic field [see Eq. ®]. 
This happens both in the case of the quasilocal anisotropy considered here and in the case where an anisotropy 
is constant along the radius (cf. Ref. [Hj]). I n contrast to this, for configurations supported by an isotropic 
fluid, the deformations due to the presence of an axisymmetric magnetic field do not affect the structure of the 
magnetic field. 

(ii) Distribution of the internal magnetic field depends considerably on the degree of fluid anisotropy and the sign 
of the parameter a. Along the radius of the configuration, the magnitudes of the components of the magnetic 
field strength Bf and Bq may be either greater or less than those of the isotropic case. In particular, in the 
case of large negative a they can differ in order of magnitude (see Fig.[l]). 

(iii) Distribution of the current is also changed considerably depending on the value of a. In particular, this is evident 
from the fact that for some values of a the current may have variable sign (see Fig. [T]l, which is impossible in 
the isotropic case. 

Therefore, we see that the structure of the internal magnetic field may depend substantially not only on the physical 
properties of matter (as demonstrated, for example, in Refs. am for neutron matter modeled by various EOS) but 
also on the degree of fluid anisotropy. 

It is clear that the obtained results are essentially model dependent and are determined by a specific manner of 
modeling the anisotropy in the system. Unfortunately, at present there is no fully reliable way to determine the true 
nature of the anisotropy and how large it may be in realistic superdense configurations. However, if high-density 
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matter of compact objects may possess an anisotropic pressure, one might expect that regardless of the way in which 
the anisotropy is modeled, its presence will result in changes of magnetic fields, which can be evaluated by using the 
approach employed here. 
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